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A MOVING-STRIP  FOURIER  ANALYZER 


by 

H.  J.Grenv  ill  e-Wells 


Laboratory  for  Insulation  Research 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts 

Abstract:  A device  embodying  Robertson's  principle  of  movable  strips  is  des- 

cribed for  performing  two-dimensional  Fourier  summations  of  the  type 

2 2 F(hk) [cos  2-rr(hX  + kY)lwithout  first  rearranging  them  in  the  form 

H K 

2 £ F(hk)  fcos  2iThX  • cos  2irkY  — sin  2irhX  • sin  2irkY  The  same 

H K 

strips  can  be  used  to  evaluate  sine  summations  and  some  three-dimen- 
sional summations.  Alternatively,  the  summations  can  be  completed 
by  using  masks.  The  construction  of  the  device,  which  can  accom- 
modate values  of  F(hk)  = +_  ( 1 — 99),  h = —15  to  15,  k = 0tol5  for  X, 

Y = l/36  of  the  repeat  distance,  is  described:  and  its  operation  illus- 
trated by  an  example. 


1 . Introduction 

This  analyzer  is  intended  to  handle  expressions  of  the  form 

H K 

S 2 F(hk)  2,(hX  + kY) 

-H  0 

directly,  without  first  rearranging  them  in  the  form 

H K , 

£ £ F(hk)[_cos  2tt1iX  * cos  2-rrkY  — sin  2whX  ’ sin2TrkYJ. 

-H  0 

A device  embodying  this  principle  was  described  by  Robertson,^  and  a 
photograph  given  showing  his  moving -strip  calculator  set  up  for  a two-dimen- 


1)  J.  M.  Robertson,  Phil.  Mag.  21,  176  (1936). 
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sional  summation.  Our  analyzer  differs  from  Robertson's  mainly  in  that  it  has 
been  constructed  in  a form  suitable  for  comparatively  inexpensive  production, 
and  in  that  special  provision  is  made  for  recording  and  using  subtotals.  Further, 
within  the  limits  imposed  by  inability  to  accommodate  more  than  16x31  = 496 
strips  at  one  time  in  a single  outfit,  a three-dimensional  summation  of  the  form 

2 2 2 F(hk£)  ^ns  2ir(hX  + kY  + fz) 

can  be  handled  in  exactly  the  same  way  as  a tv,  o-dimensional  summation. 

As  a compromise  between  resolution  and  compactness,  the  trigonometric 
functions  are  given  at  intervals  of  10°  (i.e.,  1 / 3 6 of  the  cell  edge).  In  this  case, 
the  value  at  ir/2  is  explicit,  so  that  the  same  strips  can  be  used  for  both  cosine 
and  sine  summations.  The  amplitudes  of  the  strips  range  from  —99  to  +99. 
Values  of  h = —15  to  +15,  and  of  k = 0 to  1 5 are  availab’e. 


2.  The  Basic  Principle  Employed  in  the  Moving -Strip  Method 
In  the  discussion  which  follows,  the  units  of  X,  Y and  Z are  1 / 3 6 of  the 
repeat  distance.  Since  all  possible  values  of  F(h00)  cos  2uhX  are  contained  on 
a strip  running  from  0 to  2tt,  the  first  operation  is  the  location  of  the  strips  so 
that  the  appropriate  value  is  selected  for  the  particular  point  (X)  at  which  the 
summation  is  being  evaluated.  As  an  example,  consider  the  strip  of  amplitude 
93: 


93 


93  92  87  81  71  60  46  32  16  0 IF  17  TG  . . . etc., 


where  hX  = 0 1 23456789  10  11  12  ...  etc.  (I) 

If  h = 1,  then,  as  X increases,  successive  values  of  93  cos  2irhX  will  be  those 
printed  on  the  strip;  if  h = 2,  then  successive  values  of  cos  2rrhX  will  be  for 

X=  0.1.  2.  3.  4.  5.  6...  etc.,  (II) 

dliu  w licit  h — 3 , for 

X=  0..1..2..3..4...  etc.  (Ill) 
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Fig.  1.  Basic  principle  for  the  summation  of  a one -dimensional  series  of  the 
H 

form  F(hk)  cos  2u(hX  + kY)  when  k and  X are  kept  constant.  Here  k - 1, 


X = 2. 


Hence,  if  an  arbitrary  origin  (i.e.,  M in  Fig.  1)  is  chosen,  and  a set  of 
scales  like  (I),  (II),  (III)  are  marked  off  from  this  origin,  reading  to  the  left,  the 
value  of  F cos  2nhX  can  be  obtained  as  follows:  Suppose,  as  in  Fig.  1,  that  ampli- 
tude F(3l)  = 93.  The  strip  "93"  is  chosen  and  inserted  in  the  slot  corresponding 
to  h = 3.  Suppose  that  the  electron  density  is  to  be  evaluated  at  the  point  (X,  Y) 

= (2,  5).  Since  X - 2,  the  left-hand  end  of  the  strip  is  set  against  the  number  "2" 
on  the  location  scale,  and  the  value  of  F(31)  cos  2Tr(hX  + kY)  will  appear  under 
the  number  "5"  (since  Y = 5)  on  the  right-hand  side  of  the  origin  M.  Thus  the 
value  of  93  cos(2ir/36)(3  x 2 + 1x5)  = TI.  It  will  be  seen  that  if  the  one-dimensional 
summation  having  constant  k,  such  as 

H 

£ F(hl)  cos2ir(h  • X + 1 • Y)  (IV) 

0 

is  to  be  evaluated  at  the  point  (X,  Y)  = (2,  5),  the  ends  of  the  strips  of  ampli- 
tudes F(h£)  will  lie  along  the  line  PQ,  thus  automatically  displacing  each  strip 
a distance  hX  (i.e.,  2h)  to  the  left  of  the  origin.  The  value  of  kY  at  any  point 
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(X,  Y)  is  constant  for  this  group  of  terms  of  constant  k,  and  hence  the  value  of 
one -dimensional  summations  like  (IV)  at  the  point  (2,  5)  is  obtained  by  summing 
the  numbers  in  the  vertical  column  R.  Without  altering  the  setting  of  the  strips 
(i.e.,  keeping  the  value  of  X constant)  this  one -dimensional  summation  having 
constant  k can  be  obtained  for  all  values  of  Y (i.e.,  at  all  the  points  (2,  Y)  ). 

A two-dimensional  cosine  series  can  therefore  be  handled  as  a set  of 
one -dimensional  summations,  each  of  constant  k;  and  if  these  are  arranged  on  a 
set  of  boards  each  bearing  a k location  scale  on  the  right-hand  side  (which  scale 
will  be  the  mirror  image  of  the  corresponding  h location  scale  on  the  left-hand 
side)  then  the  strips  having  amplitudes  corresponding  to  all  the  observed  values 
of  F(hkO)  can  be  inserted  on  the  boards  at  the  beginning  of  the  operation  and  used 
throughout  the  preparation  of  the  projection. 

A calculator  embodying  this  principle  must  therefore  meet  the  following 
requirements:  (l)  As  can  be  seen  from  Fig.  1,  a separate  board  will  be  needed 
for  each  value  of  the  Miller  index  k,  bearing  the  appropriate  k-location  scale 
along  the  top  edge,  while  the  set  of  h-location  scales,  which  will  be  called  the 
Location  Sheet,  will  be  the  same  for  all  boards;  (2)  Since  values  of  the  cosine 
function  must  appear  in  all  columns  like  R in  Fig.  1,  even  when  the  end  of  the 
strip  is  located  at  the  left-hand  side  of  Fig.  1,  a strip  must  carry  two  complete 
cycles  of  the  cosine  function;  (3)  Some  kind  of  guide  must  be  provided  to  keep 
the  strips  in  place;  (4)  If  the  cosine  strips  are  to  be  used  for  sine  summations, 
it  must  be  possible  to  move  the  location  sheet;  thus  two  complete  cycles  of  these 
h scales  will  be  needed,  as  for  the  strips;  (5)  Provision  must  be  made  for  re- 
cording the  totals  of  the  one-dimensional  summations,  preferably  in  a form  which 
will  be  convenient  for  future  use. 

3 . Details  of  Construction 

In  order  to  obtain  clear  and  easily  legible  numbers,  the  original  strips 
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and  location  scales  were  prepared  on  an  I.B.M.  "Executive"  electrical  type- 
writer; the  resulting  appearance  of  the  printed  version  can  be  seen  in  Fig.  1. 
Available  facilities  fixed  the  optimum  size  of  the  printed  cards  at  17  in.  x 11  in., 
which  is  therefore  the  size  of  the  location  sheet,  while  the  cosine  strips  are  17  in. 
long.  Such  a strip,  even  when  cut  from  heavy  card,  will  be  too  flimsy  if  it  is 
less  than  l/4  in.  wide  and  it  was  not  possible  to  support  for  easy  movement 
more  than  3 such  strips  per  inch.  The  working  height  of  the  location  sheet  being 
approximately  10  in.,  the  maximum  number  of  strips  was  thereby  fixed  at  31, 
i.e.,  h = —15  to  +15.  This  limitation  is  unfortunate,  since  in  fact  18  orders  should 
be  available  if  the  cell  edge  is  divided  into  36  parts. 


Fig.  2.  Layout  of  the  boards  used  for  summations. 

To  meet  requirements  (1)  to  (5),  the  design  shown  in  Fig.  2 was  adopted. 
The  backing  board  is  a 12  in.  xl8  in.  sheet  of  white-coated,  pulp  board  to  which 
the  other  components  are  affixed  with  staples.  The  co-ordinate  strip  and  the  lo- 
cation sheet  are  printed  in  red  on  high-grade  rag  card  (175  lb&).  The  Miller  in- 
dex strip  and  the  cosine  strips  are  printed  in  black  on  a heavier  card  (200  lbs.) 
of  the  same  type.  The  cover  sheet,  which  prevents  locution- sheet  numbers  ap- 
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pearing  behind  the  cosine  strips,  is  a 24-lb.  white  paper.  Requirement  (3)  is 
met  by  making  the  cosine  strips  slide  in  slots  stamped  in  2S0  aluminum  sheet, 
1/2  in.  wide  and  0.012  in.  thick,  by  means  of  a rotary  stamp  (Fig.  3).  Require- 
ment (4)  is  met  by  sliding  the  location  sheet  under  the  Miller  index  strip  and 
under  the  slotted  strip  stapled  to  the  backing  board  at  A and  B,  and  at  C and  D, 
respectively.  Requirement  (5)  is  met  by  providing  a spare  slot  for  a 1/4  in. 
strip  just  below  the  K scale. 


There  are  16  boards  (k  = 0 to  15)  of  the  type  described,  5 sets  of  strips 
having  positive  and  negative  amplitudes  1 to  33  and  3 sets  of  strips  having  + and 
— amplitudes  34  to  99,  a total  of  726  strips,  which  are  already  cut.  One  complete 
set  of  strips  of  positive  and  negative  amplitudes  1 to  99  are  laid  out  in  flat  files, 
similar  to  that  shown  in  Fig.  2 of  Robertson's  paper. ^ 


2)  J.  M.  Robertson,  J.  Sci.  Instr.  25,  28  (1948). 
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4.  Evaluation  of  a Two-Dimensional  Projection 
As  an  example  of  the  method,  the  electron  density  in  the  two-dimensional 
[010]  projection  of  a 1,  4-dimethoxybenzene  molecule  is  computed  from  data 

3) 

published  by  Goodwin,  Przybylska  and  Robertson.  ' The  space  group  is  Pbca, 

so  that  (ho£)  reflections  occur  only  when  £ = 2n.  The  observed  F(hoZ)  data 

contained  in  Table  5 of  their  paper  have  been  used  to  derive  the  values  shown 

in  Table  1,  which  are  the  amplitudes  actually  used  in  the  summation 
H L.  0 L 

D(X,  0,  Z)  = n F(hOZ)  cos  2ir(hX  + £ Z)  + £ £ F(hOi)  cos  2*(-hX  + £ Z).  In 
‘ 00  -HO 

the  ensuing  discussion  and  in  Figs.  4 and  6 only  the  first  term  of  this  expression 
has  been  used  for  illustrative  purposes,  but  Fig.  7b  and  hence  Table  2 were  ob- 
tained from  the  complete  expression. 

o o 

This  example  was  chosen  because,  in  the  unit  cell,  a = 7.29A,  c = 16.55A, 

so  that  reflections  were  observed  for  which  £ ^ 20;  hence  it  was  possible  to  see 

whether,  for  practical  purposes,  values  of  £ > 15  could  be  neglected,  as  must  be 

done  with  this  moving-strip  calculator. 

The  origin  was  taken  at  a center  of  symmetry,  and  the  projection  plotted 

for  X =(2tt/36)(0  to  9),  Z = (2tt/36)(0  to  18).  Since  Z >X  , the  one-di- 
max max’ 

mensional  summations  of  constant  / are  evaluated  rather  than  those  of  constant 

h.  One  board  is  therefore  required  for  each  of  the  8 values  of  £ = 0,  2,  14; 

strips  of  amplitudes  given  in  the  appropriate  column  of  Table  1 are  inserted  in 
the  correct  h slots  on  each  board,  as  can  be  seen  in  Fig.  4 on  the  board  for  which 
£ - 14.  To  calculate  the  electron  density  at  a row  of  points,  say  (3,  0,  Z),  the 
left-hand  end  of  each  strip  is  placed  against  the  marker  for  "3"  on  its  h-location 
scale,  and  the  summations  performed  for  each  column  shown  to  be  necessary  by 
the  markers  on  the  £ scale.  Since  cos  (A  + ir)  = —cos  A,  it  is  never  necessary  to 
compute  more  than  18  out  of  the  possible  36  subtotals  on  any  board,  and  often 

A 

fewer  are  needed;  for  instance,  ir.  the  example  chosen,  in  which  Z = 2n,  it  can 

be  seen  from  the  scales  for  all  the  8 boards  visible  in  Fig.  4 that  only 

3)  T.H.  Goodwin,  M.  Przybylska,  and  J.  M.  Robertson,  Acta  Cryst.  3,  279  (1950). 
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Fig.  4.  (a)  Evaluation  of  the  one-dimensional  series  ap  F(h,  14)  cos  2n(h,  3 + 14. Z) 
at  the  points  (3,  0,  Z).  On  the  board  l - 14,  strips  of  amplitudes  F(h,  0,  14)  ob- 
tained from  Table  I,  are  placed  in  the  b slots.  Since  X = 3,  the  left-hand  end  of 
each  strip  is  set  against  the  number  "3"  on  the  Location  Sheet,  and  the  subtotals 
for  each  value  of  Z written  on  the  subtotal  strip  in  the  top  slot.  Such  a strip 
for  X = 3 is  prepared  for  each  value  of  /,  i.e.,  for  each  board,  (b)  Completion 
of  the  p (3,  0,  Z)  summation.  If  the  boards  are  arranged  in  a staggered  pile  as 
shown,  the  final  summation  at  any  point  (3,  0,  Z)  can  be  picked  out.  In  the  figure 
the  indicators  are  set  to  pick  out  the  summation  at  (3,  0,  2). 


i 


T.R. 64 


-10- 


l+9+9  + 3+  9 + 9 + 3 + 9=  52  subtotals 
have  to  be  computed  for  each  value  of  X / 0.  The  strips  must  be  reset  each  time 
X is  changed,  i.e.,  for  every  line  of  constant  X in  the  projection. 

5 . Methods  of  Completing  the  Two-Dimensional  Summation 

A.  Evaluation  of  /P(X,  0,  Z)  at  each  point  in  turn.  One  especially  useful 
feature  of  this  method  is  that  the  electron  density  at  any  selected  point  (X,  0,  Z) 
can  be  obtained  independently  by  setting  all  the  strips  for  the  value  "X"  on  the 

h scales,  and  adding  together  the  "Z"  column  from  each  of  the  8 boards  without 
prior  evaluation  of  subtotals.  If  the  electron  density  is  required  at  a small 
number  of  arbitrarily  chosen  points  in  the  projection,  this  method  is  very  fast. 

If,  however,  the  values  of  yb(X,  0,  Z)  are  required  for  the  complete  projection, 
it  is  better  to  make  use  of  subtotals  for  two  reasons: 

1.  For  every  value  of  l which  is  a factor  of  36,  or  contains  a 
factor  of  36,  the  same  column  of  numbers  will  be  used  for  two  or  more  values 
of  Z.  For  example,  when  S - 12,  the  three  columns  indicated  in  Fig.  4 are  each 
used  for  12  values  of  Z.  It  is  a considerable  waste  of  time  to  sum  these  columns 
again  for  each  of  these  12  values,  especially  whe  i there  are  a large  number  of 
strips  on  the  board. 

2.  An  appreciable  time  is  inevitably  consumed  in  the  manipulation 
of  strips  and  boards,  and  this  is  reduced  to  a minimum  by  using  subtotals. 

B.  Use  of  subtotal  sheets.  In  Fig.  4 it  can  be  seen  that  when  the  sub- 
totals have  been  computed  for  a given  value  of  X,  the  final  summation  at  any 
point  (X,  0,  Z)  will  be  obtained  by  yimming  the  appropriate  subtotals,  which 
can  be  picked  out  on  the  subtotal  strips.  In  Fig.  4 the  subtotals  required  for 

Z = 2 are  shown  by  ring  indicators.  It  is  quite  possible  to  complete  a summa- 
tion in  this  way,  but  the  staggered  pile  of  boards  becomes  unwieldy  when  ail  16 
boards  are  in  use;  hence  a subtotal  sheet  (reproduced  in  Fig.  5)  containing  these 
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Fig.  5.  Location  scales  for  0,  1, 2.  . . . 15.  When  all  16  boards  are  in  use,  the  stag- 
gered pile  (Fig.  2)  becomes  unwieldy.  If  a sheet  of  tracing  paper  is  superposed  on 
Fig.  5.,  the  subtotals  for  a particular  value  of  X can  be  entered  on  it,  and  the  final 
totals  prepared  from  this  instead. 
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scales  is  supplied  with  the  device:  the  subtotals  can  be  written  on  a sheet  of 
tracing  paper  superposed  on  these  scales,  and  the  final  totals  prepared  from 
this  instead;  the  tracing  can  then  be  kept  for  reference. 

Keeping  such  a record  of  the  subtotals  for  each  value  of  X has  an  ad- 
vantage if,  in  successive  Fourier  summations,  the  sign  of  an  F(hO£)  term 
has  to  be  changed,  and  this  is  done  by  adding  a term  of  opposite  sign  having 
tv/ice  the  amplitude  of  F(ho£).  The  appropriate  cosine  strip  can  be  inset, 
in  the  second  spare  slot,  directly  beneath  the  subtotal  strip,  with  its  left-hand 
end  located  at  the  required  value  of  X in  the  usual  way,  and  then  the  corrected 
subtotals  are  obtained  merely  by  adding  the  two  strips  together. 

This  designation  of  subtotals  by  means  of  ring  indicators,  as  in  Fig.  4, 

2 5) 

is  equivalent  to  the  use  of  masks  for  the  selection  of  terms.  ’ ' Although 

they  could  easily  be  made,  no  masks  are  provided  with  this  device  because  it 
is  felt  that  while  they  would  certainly  reduce  the  probability  of  error,  masks 
are  not  absolutely  necessary  as  the  value  of  Z is  explicitly  stated  above  each 
subtotal. 

C.  Subtotal  strips  treated  like  cosine  strips.  On  the  board  for  which 
|=0,  there  will  be  space  for  strips  of  amplitudes  F(h00)  and  F(K00).  Since,  in 
general,  F(h00)  = F(K00),  only  one  of  these  groups  will  be  needed,  and  either  set 
of  slots  may  be  used. 

If,  therefore,  the  F(K00)  slots  are  used,  as  in  Fig.  6(a),  then  the  F(h00) 

slots  will  be  available  for  use  with  the  subtotal  strips.  All  the  subtotal  si  ’ps 

visible  at  the  top  of  Fig.  4 arc  placed  on  the  board  for  which  1 = 0,  the  strip 

from  the  board  t - n being  placed  in  the  slot  for  which  h = n.  The  summation 

at  all  the  points  (3,  0,  Z)  in  turn  is  obtained  by  setting  the  strips  against  the 

number  "Z"  on  the  location  sheet,  and  summing  the  numbers  in  the  first  column 

to  the  right  of  the  origin  (M  in  Fig.  1).  The  strips  in  Fig.  6 are  set  to  give  the 

°2 

value  of  p{ 3,  0,  2)  = 807,  which,  as  stated  in  Table  1,  is  equivalent  to  8.07e/A  . 
5)  A.  I_.  Patter  son  and  G.  Tunell,  Amer.  Min.  2_ 7,  655  (1942). 
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Because  a certain  amount  of  time  is  inevitably  consumed  in  manipulating 
the  boards  and  strips,  it  is  difficult  to  estimate  the  time  taken  for  a summation. 
The  relative  speeds  of  the  alternatives  already  described  depend  partly  on  the 
number  of  strips  in  use  on  each  board,  and  partly  on  the  shape  of  the  section  of 
the  total  projection  which  has  to  be  computed.  For  example,  when  subtotals 
are  used,  if  all  the  possible  subtotals  for  a fixed  value  of  X are  computed  for 
each  board  in  turn  (the  fastest  method  of  subtotalling),  this  will  in  general  have 
to  be  done  regardless  of  whether  Z runs  from  0 to  35,  or  only  l/2  or  l/4  this 
range.  Also,  it  has  to  be  done  again  for  each  value  of  X,  so  that  the  subtotalling 
method  is  much  faster  (indeed,  nearly  twice  as  fast  if  the  boards  are  well- 
filled)  when  the  projection  has  to  be  evaluated  at,  say,  Z = 0—35,  X = 0—8,  than 
when  it  must  be  done  for  Z = 0—17,  X = 0—17,  although  the  summation  is  required 
at  the  same  number  of  points. 

6.  Accuracy  and  Applicability  of  the  Method 

The  completed  projection  obtained  from  this  summation  is  shown  in  Fig. 

3) 

7,  together  with  the  projection  published  by  Goodwin  et  al.  ' for  comparison. 

The  general  shape  of  the  molecule  is  easily  recognizable,  and  the  co-ordinates 

of  the  atomic  centers,  taken  merely  as  the  estimated  positions  of  the  highest 

o 

points  of  the  peaks  on  a contour  map  on  the  scale  1 inch  = 1A,  are  given,  to- 
gether with  the  published  values,  in  Table  2.  In  view  of  the  fact  that  some  of 
the  high-order  terms  neglected  were  quite  large  (up  to  13  percent  of  ^max)»  an<* 
that  no  attempt  was  made  to  estimate  the  true  centroid  of  a peak,  the  excellent 
agreement  obtained  in  some  cases  is  probably  fortuitous. 

This  raises  the  following  question:  suppose  atomic  co-ordinates  esti- 
mated from  a Fourier  projection  on  a rather  coarse  mesh  (as  has  been  done 
here)  without  elaborate  attempts  to  obtain  high  accuracy,  are  used  to  compute 
structure  factors,  and  that  these  give  a value  of  R ^ 30Jl ; will  this  be  a satis- 
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factory  basis  for  the  application  of  a least- 
squares  refinement  process?  Work  now  in 
progress  on  the  least-squares  method  sug- 
gests that  this  will  be  true  in  at  least  a cer- 
tain number  of  cases.  If  so,  a Fourier  summa- 
tion device  of  the  type  described,  with  values 
obtainable  at  10°  intervals,  would  have  a wider 


application  than  might  otherwise  be  expected. 

In  any  case  it  will  be  adequate  for  work  on  crys- 
tals in  which  reflections  with  h,  k,  t >15  are 
not  observed. 


Fig.  7.  [010]  projection  of  1,4- 
dimethoxybenzene  molecule:  (a) 
from  data  of  Goodwin  et  al.  ' ; 
(b)  as  computed  from  data  in 
Table  1,  where  l $ 14. 


7.  Evaluation  of  Series  Other  than  the 
Two-Dimensional  Cosine  Series  In- 
volving cos  2tr(hX  + k.Y). 

Two  important  extensions  of  the  method 


are  indicated  in  the  caption  to  Fig.  8.  In  addition  to  these,  summations  of  the 
form  cos  2it(hX  — kY)  can  be  obtained  by  placing  the  F(hk)  strips  in  the  h slots, 
since  cos  2ir(— hX  + k Y)  = cos  2ir(hX  — kY). 

Another  possibility  is  the  inclusion  of  phase  angles.  Robertson^  sug- 
gested that  for  noncentro symmetric  projections,  am  additional  displacement  of 
the  strip  corresponding  to  the  phase  amgle  could  be  used,  with  the  limitation  that 
it  can  be  varied  only  in  discrete  amounts  (multiples  of  10°  in  this  analyzer). 

It  is  also  possible  to  carry  out  the  Beevers-Lipson  type  of  summation^ 
involving  expressions  like  F(hk  t)  S*n  2trhX  * S*n  2irkY  • S*n  2irlZ  by  evaluating 
F(hk  l)  2irhX,  and  using  this  as  a new  aunplitude  F'(hk/),  and  so  on.  In  this 
case  the  one -dimensional  summations  required  can  be  obtained  by  locating  the 
strips  as  usual  and  evaluating  cos  2ir(hX  + kY),  where  Y = 0.  Alternatively,  if 


6)  C.  A.  Beevers  and  H.  Lipson,  Proc.  Phys.  Soc.  48,  772  (1936). 
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Table  2.  Co-ordinates  x and  z of  atoms  in  1,  4-dimethoxybenzene  molecule. 


Atom 

G.  P.&  R. 

X 

El 

G.  P.&R. 
z 

G-W 

z' 

HH 

C(l) 

0.0055 

0.  008 

-0. 002 

0.0777 

0.  082 

-0. 004 

C(  2) 

0.0927 

0.  074 

+0. 016 

0.0612 

0.  061 

0 

C(  3) 

0.0977 

0.10« 

-0.  011 

-0.0156 

-0.  021 

+0. 005 

0 

0.0104 

0.  010 

0 

0.1541 

0.164 

-0.  010 

C(CH3) 

-0.0763 

-0. 069 

+0. 007 

0.1730 

0. 166 

+0. 007 

all  the  strips  were  set  for  X * 0,  they  would  present  the  same  appearance  as 
Robertson's  sorting  board, except  that  values  of  the  cosine  function  are  listed 
for  the  whole  cycle  0 to  Ztr  instead  of  only  0 to  tt/2,  so  that  the  sign  of  each  term 
is  explicit  and  does  not  have  to  be  supplied  by  auxiliary  markers.  It  would  there- 
fore be  possible  to  pick  out  the  summation  by  placing  the  strips  in  the  subtotal 
slots,  and  proceeding  as  in  section  5(b),  or  by  making  masks.  However,  since 
devices  are  obtainable  which  have  been  specifically  designed  for  computing 
Fourier  summations  arranged  as  sine  and  cosine  products,  this  represents  a 
need  which  has  already  been  satisfactorily  met,  and  therefore  no  masks  have 
been  provided. 
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